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Abstract 

We show on the example of the discrete heat equation that for any given discrete 
derivative we can construct a nontrivial Leibniz rule suitable to find the symmetries 
of discrete equations. In this way we obtain a symmetry Lie algebra, defined in 
terms of shift operators, isomorphic to that of the continuous heat equation. 



1 Introduction 



Lie point symmetries were introduced by Sophus Lie for solving differential equations. 
They turn out to provide one of the most efficient methods for obtaining exact analytical 
solutions of partial differential equations This essentially continuous method has 

been recently extended to the case of discrete equations 

Let us write a general difference equation, involving, for notational simplicity, one 
scalar function u{x) of p independent variables ) evaluated at a finite 

number of points on a lattice. Symbolically we write 

E{x, T'^uix), T'''A,Mx), r'=-A,^A,^,u(a;), . . .) = 0, (1.1) 

where E is some given function of its arguments, 

TMa^):={^;;T;^■■T;;^.(x)^ a = {a,,a,, . . . ,a,), z = l,2,...,p, (1.2) 

with ai,mi,ni fixed integers (mj <ni), and 

T'^'uiyx) = u{xi,X2, . . . , Xi + ttiai, x^+i, . . . , Xp), (1.3) 

(the other shift operators T^', T'^^^ are defined in a similar way), A^^- is a difference 
operator which in the continuous limit goes into the derivative and ai is the positive 
lattice spacing in the uniform lattice of the variable Xj (i = 1, . . . ,p). 

The simplest extension of the continuous case is when the symmetry transformation 
which leaves the equation on the lattice invariant depends just on u{x) and not on its 
shifted values, what was called intrinsic point transformations [§]. In this case the whole 
theory goes through but the resulting transformations are somehow trivial and provide 
often not very interesting solutions. 

In Ref. [0 it has been proved that the intrinsic transformations can be extended for 
linear equations by considering symmetries where the transformed function u{x) depends 
not only on the old u{x) and x but also on the function u in shifted points T°'^u{x) on 
the lattice. 

The results presented in |^ where obtained considering the difference operator 

A.. = A+ = (1.4) 

This difference operator is the simplest one which, when Uj — 0, goes into the standard 
right derivative with respect Xi. 

Obviously, more general definitions of the difference operator can be introduced and 
one would like to be able to prove that the obtained symmetries are independent on the 
discretisation of the difference operator one is considering. Among the possibilities, let us 
mention 

1 - T-i 



corresponding to the left derivative, and the symmetric derivative (which goes into the 
derivative with respect to Xi up to terms of order af) 

rp _ rp-l 

K, = -^^^ (1.6) 

Using the techniques of approximate differentiation ||^ one could write more complicate 
formulas for difference operators which, however, are out of the scope of this letter. 

In the following we will show that by a nontrivial definition of the Leibniz rule we 
can construct symmetries for any difference operator. Section ^ is devoted to review the 
Lie group formalism for difference equations. In Section |^ we apply it to the case of the 
discrete heat equation using the discrete derivatives ( p..5| )-( pr6D , finding difficulties in the 
case ( p..6|) . In Section ^ we show via examples how one can introduce in a consistent 
way the Leibniz rule and use this result to obtain determining equations equivalent to 
those for the standard discrete derivatives (i.e., A^.). This fact shows that we get differ- 
ent representations of the same symmetry group. In the Conclusions we sketch how to 
write the determining equations when one uses a generic difference operator by using the 
corresponding Leibniz rule obtained with our approach. 



2 Lie symmetries of difference equations 

Among the different algebraic methods for calculating the symmetries of discrete equations 
we will make use in this paper of the approach presented in Ref. [|l^, based on 
the formalism of evolutionary vector fields [Q. 

For a difference equation of order like that given in (|1 . 1|) the infinitesimal symmetry 
vectors in evolutionary form, which in the continuous limit go over to point symmetries, 
take the general expression 

X, = Qdu = ^,(x, T'^u, a,, at)T'A^M - 0(x, T'^u, a,,, at)j du, (2.1) 

with ^i{x, T'^u, ax, crt) and (j){x, T^u, a^, at) functions which in the continuous limit go over 
to ^i{x,u) and (f){x,u), respectively. 

The vector fields Xe generate the symmetry group of the discrete equation (|1.1D, 
whose elements transform solutions u{x) of the equation into solutions u{x). The iV-th 
prolongation of X^ must verify the invariance condition 

pr^XeE\E=o = 0. (2.2) 

The group generated by the prolongations also transforms solutions into solutions, and 
A^^u, A^^A^ u, . . . (up to order N) into the variations of u with respect to Xj. The formula 
of pr'^Xg is given by 

pr^Xe = J2T''QdT^u + E^^^Q'^^T^^A.^^ + E^''^Q'''^^T^-A.,A... + ■ ■ ■ (2-3) 
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The summations in ( |2.3| ) are over all the sites present in ( [I.ID , Q^^^^ , . . . are total 
variations of Q, i.e., 

Q-"=Aj^g, Q^^^^=AlAlQ, (2.4) 

where the partial variation A^- is defined by 

A^f{x,u{x),A^u{x),...) = ^[f{x + a,u{x),{A^u){x),...) 

-f{x,u{x),A^u{x), . . .)], x = Xi, 

and the total variation A^. by 

Alf{x,u{x),A^u{x),. . .) = ^[f{x + a,u{x + a),{A,^u){x + a),...) 

(2.6) 

-f{x,u{x),A^u{x),...)], Wx = Xi. 

The symmetries of the equation ( |1 . 1|) are given by equation ( p.2|) . From it we get the 
determining equations for and as the coefficients of linearly independent expressions 
in the discrete derivatives T^A^-^m, T^^A^.^m, .... 

The Lie commutators of the vector fields are obtained by commuting their first 
prolongations and projecting onto the symmetry algebra Q, i.e., 



dQ2 /O dQi.^XidQ2 r\^i dQi \ a (^■'^) 

where the du^^ terms disappear after projection onto Q. 

The formalism presented above is very complicate and if the system is nonlinear is 
almost impossible to get a result since the number of terms to consider is a priori infinite. 
The situation is simpler in the case of linear equations for which we can use a reduced 
Ansatz. We assume that the evolutionary vectors ( p.l| ) have the form 

(2.8) 



Now the vector fields X^ can be written as Xe = {Xu)du, i.e., 

^ = E6(2:,T^a,,a^)A,, -0(x,r",a,,at). (2.9) 

i 

The operators X may span a subalgebra of the Lie symmetry algebra (see Ref. fTD])- 

3 Discrete heat equation 

Let us consider the equation 

{At - A,,)u{x) = 0, (3.1) 
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which is a discretisation of the heat equation. As the equation is hnear we can consider 
an evolutionary vector field of the form 

X, = Qdu = {rAt + CA,u + fu)du, (3.2) 

with T, ^ and / arbitrary functions of x,t,Tx Tt, ax and at- As and Tt are operators 
not commuting with x and t, respectively, r, ^ and / are operator valued functions. 



Since equation (|3.1|) is a second order difference equation it is necessary to use the second 



prolongation. The determining equation is 

AjQ - ALQ|a..«=a.. = 0, (3.3) 

which explicitly reads 

At(CA,M)+At(rAtM)+At(/M)-[A,,(eA,M)+A,,(rAtM)+A,,(/M)]|A,,«=A,« = 0. (3.4) 

We had no need to give the explicit form of the operator A to get equation ( p.4| ). Only 
when developing expression ( p.4|) we need to apply Leibniz rule and, hence, the results 
will depend from the given definition of the discrete derivative. 

3.1 Symmetries in the right (left) discrete derivative case 

Choosing as in Ref. 0, [1^ the derivative A"*" and consequently the Leibniz rule 

A^{fg)=A^{f)Tg + fA+g (3.5) 
we obtain the following set of determining equations 

A+r = 0, 

(A+r)T, - 2{AtOT. = 0, 

{AtoTt - {AU)n - mtm = o, 
{Atm - {^tj)n = 0, 

by equating to zero the coefficients of A^tU, AtU, A^u and u, respectively. The solution 



of (|3.6|) gives 

T = t(2)r2 + tTi + To, 

e = |x(ri + 2tr2)TtT-^ + t^^ + ^o, (3.7) 

where tq, ti, T2, ^0, ^1 and 7 are arbitrary functions of T^, Tt and of the spacings ax and 
at, and x^"'\ t*-"-* are the Pochhammer symbols given by, for instance, 

= x{x - a^) . . . (x - (n - l)a^). (3.8) 
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a 


suitable choice of the functions Tj, ^i, 


and 7 we 


get the followin 


Po 


= {Atu)du, 




(ro = 1) 


Pi 


= {A,^u)du, 




(^0 = 1) 


W 


= udu, 




(7 = 1) 


B 


= {2tTf^A,u + xT-^u)du, 




(6 = 2rr^) 


D 


= {2tTt-^Atu + xT-^A^u + lu)du, 




(ri = 2rr\7 


K 


= {t^Tf'^Atu - attTf^Atu + txTf^T^ 


-^A^u 






+ \x^T-^u - \a.,xT-^u + itrr 




(^2 = T,-') 



(3.9) 



that close into a 6-dimensional Lie algebra which is isomorphic to the symmetry algebra 
of the continuous heat equation. 

A second choice for the discrete derivative is A~ and now the Leibniz rule becomes 

A~U9) = A-U)T-'g + fA-g (3.10) 
This gives the same results ( p.7|) and ( p.9|) provided we make the substitution T — > T^^. 

3.2 Symmetries in the symmetric discrete derivative case 

Next let us consider the case of the symmetric derivative A*. It seems natural to propose 
the following Leibniz rule 

A^(/^) = A\f)T-^g + {Tf)A'g. (3.11) 

In this case taking as before coefficients of the different discrete derivatives of u {Auu, 
AxtU, AtU, AxU and u, respectively) the determining equations are 

T,T - T^T = 0, 

r,e-r2^-2(r.A^.r)T-i = o, 

(A?r)Tr^ + (T,/) - 2{TxA%0T-' - (ALr)T-2 - T^f = 0, (3.12) 
{AtOTf' - {AIJ)T-' - 2(r.A^/)T-i = 0, 
{Atm-' - {AU)T-' = 0. 

Note that there is one equation (the first one, associated to Auu) more that in the two 
previous cases. This implies that the solution of such equations ( ^.12[ ) is just 



r = To, ^ = ^0, / = /o, (3.13) 
with To, and /o are arbitrary functions of T^, Tt and of the spacings and ctj. 

Obviously, in this last case our naive approach does not allows us to recover the whole 
symmetry algebra of the heat equation ( [3.91 ). In fact, we will show in next Section that 
the decomposition of equation ( |3.4| ) into equations (|3.12|) via Leibniz rule (|3.11|) is not 
the most appropriate one. 
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4 Leibniz rule and symmetries with symmetric 
derivatives 



Let us consider a formal way to get Leibniz rule in the continuous case which is easily ex- 
tendible to the case of discrete derivatives. This result will allow us to get the symmetries 
of a discrete equation independently of the discrete derivative used. 

Starting from the well known result [S^i^;] = 1 by algebraic methods we get that 
[dx, /(a;)] = .f\x) (at least for analytic functions). Consequently, the Leibniz rule d^ifg) = 
fg' + f'g can be derived. 

Now, let us consider the commutator with the differential operator substituted by a 
difference one. By direct computation we have 



[Alx] = T^^\ [A:.,x] = ^^^. (4.1) 

By introducing a function (5x = (3{Tx) we can always rewrite the commutation relations 
(EID as 

[A,,x/5J = L (4.2) 
Let us note that A^Px = PxAx and At(3x = PxAt- 

For the standard left and right derivatives A^^ we easily find that = T^^, and we 
get that 

[A±,/(x)/?±] = (A±/(x)). (4.3) 
From this last expression we recover Leibniz rules (|3.5|) and ( p.lO|) . 

For the symmetric discrete derivative A^ we obtain the formal expression 

f3: = 2iTx + T-Y\ (4.4) 

and 

[A:, f{x)P:] = (Alfix)) TxP: + (T^T'fix) - fix)) PlAl. (4.5) 
From relation (|4.5| ) we get the Leibniz rule 

A^. (/(x)^(x)) = f{x)Alg{x) + [ j- ((T-i - l)/(x)) (T. - m'') 

+ {Alf{x))Tx]g{x). 

Let us note that expression ( |4.(:j ) looks very different from ( |3.11D . 



(4.6) 
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Formula ( |4.(j| ) allows us to write down the determining equations ( ^.41 ) as 
((l-T-»T-i + ((T,-l)r)T. = 0, 
l-Tf')T)Ti' + {{Tt-l)r)Z 



1 

2crt 



1 
2tTt 



1 

2(Tt 



[((1 - 7^.'^ + ((T. - 1)0 = 0, 

((1 - Tr^)e) Tf' + m - 1)0 T,] - ^ [((1 - T-1)/) T,7i + ((T. - 1)/) (4.7) 
-4^ [((1 - T-rO + 2 ((T. + T-i - 2)0 + ((T. - 1)^0 7;^] = 0, 
l-Ti')f)Ti' + iiT,-l)f)Tt 

h [((1 - T-'?f) + 2 ((T. + T-i - 2)/) + ((T,. - 1) V) T'.] = 0, 



obtained as coefficients of A^tu, AfU, A^u and u, respectively. So, we get the same number 
of equations as in the cases of the standard discrete derivatives A='=. 

The solution of eqs.( ^.7| ) is given by 

r^ = t(2)r2 + tri + ro, 

e = Ix {2tr2 + n + a,T,-'(3tT,) )-i/?|. + t^i + ^o, (4-8) 

where r2, ti, tq, ^i, and /o are arbitrary functions of T^, Tf, and aj. 

From (|4.8|) and (|3l^ ) we obtain, with a suitable choice of T2, ri, Tq, ^1, ^0 and /o, the 
following symmetries 





= (A?«)a„, 


(to = 


1) 


P( 


= (A»a„, 


(eo = 


1) 




= udu, 


(/o = 


1) 


B' 


= {2tpiAlu + xf3lu)du, 


(6 = 


m) 


D' 


= {2t(3lAlu + x^l^Alu + i^/)a«, 


(rl = 






= ((t2(/5;)2 _ tat\[3tfAl)Alu + te/?,^/?|A> 


(^2 = 






-ixCT.2(/3|)3^.„+ i^2^^s)2„+ l^^s^)^^^ 


Tl = 


= atPf-cj^,{3fAl 



(4.9) 



They close the same 6-dimensional Lie algebra generated by the operators (|3.9| ). The 
above result deserves some comments. First, there appear functions [31 {(3^) of Tt (T^) 
that can only be understood as infinite series developments. Therefore, some of the above 
symmetries ( |4.9| ) have not a local character in the sense that they are not polynomials 
in T^^,T^^. Although such symmetries give rise to the classical symmetries in the limit 
ax — >■ 0, CTf — i> 0, one of them (K) also includes surprisingly a term in (At)^ (which vanishes 
in the continuous limit since it has a factor cr|). 
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5 Conclusions 



Analysing the different Leibniz rules used in Sections ^ and we see that when we get 
the correct result the Leibniz rule must have the form 

A, {f{x)g{x)) = f{x)AMx) + D,{f{x))g{x), (5.1) 

where is a function of T^, (3x and (we should have written Dx{f{x);Tx, f^x,o'x), 
with /3a; given by eq. (|4.2|) , but for the sake of shortness we will simply write Dx{f), and 
similarly for Dt{f)). Once we have chosen a particular discrete derivative A^,. we can write 
the explicit expression of D^. In particular, for A^ and A^. their corresponding D.j.{f) 
functions are 

= ((t-i-i)/)a:. + ((a:,)/)tj 



Using the general Leibniz rule (|5.1| ) for an arbitrary discrete derivative we obtain from 
equating to zero the coefficients of A^tu, AfU, Aj.u and u, respectively, the following 
set of determining equations 

Dxir) = 0, 
A(r)-2D.(0=0, 
Dt{0-Dxx{0-2Dx{f) = 0, 
A(/) - Dxxif) = 0. 

where D,.,.(/) = Dx{Dx{f)). 

For the cases A^ and A* from (|5.3| ) we recover the determining equations ( p.6|) and 
( [4.7|) , respectively. 

We have always four equations, whose general solution will depend on the formal 
expression of the corresponding A. In fact, as long as we are able to find for the discrete 
derivative A the operator f3 satisfying the commutator (|4.2| ) this will allow us to give 
particular solutions to the determining equations ( |5.3| ) keeping the same Lie structure of 
the symmetries of the classical heat equation. Of course, some of these discrete symmetries 
can have a rather complicated expression with a nonlocal character, but this feature is a 
consequence of the structure of the operator (3. 

This procedure can be straightforwardly applied to other discretisations such as the 
wave equation [|lT| or even equations including a potential term. 



Something similar happens with discrete equations on a g-lattice: once we have a 
solution for f3 we can derive symmetry operators with a Lie structure (a situation which 
is seldom considered). On the other hand, if we are more interested in local symmetries, 
with some restrictions, the natural structure is that of a g-algebra. However, these issues 
are out of our present scope and they will be published elsewhere. 
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